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Helicity is a quadratic invariant of the Euler equation in three dimensions. As the energy, 
when present helicity cascades to smaller scales where it dissipates. However, the role 
played by helicity in the energy cascade is still unclear. In non-helical flows, the velocity 
and the vorticity tend to align locally creating patches with opposite signs of helicity. 
Also in helical flows helicity changes sign rapidly in space. Not being a positive definite 
quantity, global studies considering its spectral scaling in the incrtial range arc incon- 
clusive, except for cases where one sign of helicity is dominant. We use the cancellation 
exponent to characterize the scaling laws followed by helicity fluctuations in numerical 
simulations of helical and non-helical turbulent flows, with different forcing functions and 
spanning a range of Reynolds numbers from « 670 to 6200. The exponent is a measure 
of sign-singularity and can be related to the fractal dimension as well as to the first order 
helicity scaling exponent. The results are consistent with the geometry of helical struc- 
tures being filamentary. Further analysis indicates that statistical properties of helicity 
fluctuations do not depend on the global helicity of the flow. 



1. Introduction 

Theories of fully developed turbulence rely on the energy conservation in inviscid flows 
to derive the scaling laws of a turbulent system. However, energy is not the only quadratic 
invariant of the three dimensional (3D) incompressible Euler equation, there is also the 
helicity 



H = 



J d 3 2iv(x) • W(x), 



(1.1) 



where v(x, t) is the velocity field and w(x, t) = V x v is the vorticity field ( Moffatt |1969 1 . 
Helical motions are observed in a wide variety of geophysical flows, such as supercell con- 
vective storms which seem to owe their long life and stability to the helical nature of their 
circulation ( Lilly |[l986 1 . In astrophysics it is well known that helicity plays a key role in 



the generation of large scale magnetic fields through dynamo action ( Pouquet et al. 1976 



Moffatt 19781. In hydrodynamic isotropic and homogeneous turbulence, recent studies 
of flows with net helicity confirmed that there is a joint cascade of energy an helicity 
to smaller scales ( Borue fc Orszag ||1997 Chen et al. ||2003« Gomez &i Mininni~||2004 1 , 



as expected from theoretical arguments 



Brissaud et al. 1973 



In the case of non-helical 



flows, velocity and vorticity tend to align locally creating patches of strong positive and 
negative helicities. In helical flows, helicity still fluctuates strongly around its mean value, 
both in space and time, creating locallized patches with helicity much larger or smaller 
than the mean. As helicity is not positive definite, the study of its turbulent fluctuations 
is difficult and has received less attention. The scaling laws followed by helical fluctu- 
ations are unclear, and few comparisons are available between helical and non-helical 



2 



P. Rodriguez Imazio and P.D. Mininni 



flows to understand their impact on the energy cascade. It has been proposed that heli- 
cal structures may have an impact on the dynamics of turbulent flows as helical regions 
have the velocity and vorticity aligned, with a resulting quenching of the non-linear term 
in the Navier-Stokes equation ( Moffatt fc Tsinober 1 1992 1. Intermittency in the helicity 
flux has also been recently studied in numerical simulations ( Chen et aZ.|20036 l, and the 



geometry of helical structures was considered using wavelet decompositions (Farge et al. 
2001]) and Minkowski functionals in the magnetohydrodynamic case ( Wilkin et aZ.|200"7 1 . 

In this work we use the the cancellation exponent ( Ott et aT~||1992 l to study the fast 
fluctuations of helicity in turbulent helical and non-helical flows. The cancellation expo- 
nent was introduced to study fast changes in sign of fields on arbitrarily small scales, and 
is based on the study of the inertial range scaling of a generalized partition function built 
on a signed measure. Under some conditions, it was shown to be related to the fractal 
dimension of structures in the flow and to the Holder exponent (Vainshtein et al. 19941. It 
was used to characterize fluctuations in hydrodynamic turbulence and magnetohydrody- 
namic dynamos ( Ott et al. ||1992 1, as well as in two dimensional magnetohydrodynamic 



turbulence ( Sorriso-Valvo et al. 2002 Pietarila Graham et al. 2005). In solar wind 



observations it was used to show that magnetic helicity is sign singular (Bruno & Car 
bone |[l997 l. However, to the best of our knowledge, no studies of helical fluctuations in 



hydrodynamic turbulence have been done using the cancellation exponent. We analyze 
data stemming from direct numerical simulations (DNSs) of the Navier-Stokes equation 
in a three dimensional periodic box at large resolutions (up to 1024 3 grid points) with 
different external mechanical forcings. Both helical and non-helical flows are considered. 
We find that helicity fluctuations are sign singular in both cases, and the scaling of he- 
licity fluctuations seems to be independent of the net helicity of the flow. Furthermore, 
we obtain a relation between the helicity cancellation exponent, the fractal dimension of 
helical structures and the first order scaling exponent of helicity. 



2. The cancellation exponent 

The cancellation exponent ( Ott et al. |1992 ) was introduced to characterize the behav- 



ior of measures that take both positive and negative values, and to quantify a form of 
singularity where changes in sign occur on arbitrarily small scales. We can introduce the 
cancellation exponent for the helicity considering a hierarchy Qi(l) of disjoint subsets of 
size I covering the entire domain occupied by the fluid Q(L) of size L. For each scale I, 
a signed measure of helicity is introduced as 

Mi(0= / d 3 xH(x)/f d 3 x\H(x)\, (2.1) 

JQ t {l) I JQ{L) 

where i/(x) = v(x) • w(x) is the helicity density such that the total helicity is H 



J d 3 xH(x). Since the normalization factor in Eq. (2.1l is the integral of the absolute 
value of H over the entire domain, the signed measure is bounded betwen 1 and —1 
and can be interpreted as the difference between two probability measures, one for the 
positive component and another for the negative component of the helicity density. To 
study cancellations at a given length scale, we build the partition function 



(2.2) 
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Run 


N 


F 




pn 


V 


Re 


Rx 


Tl 


256 


TG 


2 


0.03 


2 x 10~ 3 


675 


300 


T2 


512 


TG 


2 


0.01 


1.5 x 10~ 3 


875 


350 


T3 


1024 


TG 


2 


0.03 


3 x 10~ 4 


3950 


800 


Al 


256 


ABC 


3 


0.79 


2 x 10~ 3 


820 


360 


A2 


512 


ABC 


3 


0.81 


2.5 x 10" 4 


2520 


670 


A3 


1024 


ABC 


3 


0.83 


2.5 x 10~ 4 


6200 


1100 



Table 1. Runs used for the analysis. N is the linear grid resolution, F is the forcing, either 
Taylor-Green (TG) or Arn'old-Beltrami-Childress (ABC), kF is the forcing wave number, pn 
is the relative helicity, v is the kinematic viscosity, R e is the Reynolds number, and R\ is the 
Reynolds number based on the Taylor scale. 



In the inertial range, the scaling law followed by the cancellations can be studied by 
fitting 

x(0~r«, (2.3) 

where k is the cancellation exponent. This exponent represents a quantitative measure 
of the cancellation efficiency. If such a scaling law exists and its range increases as the 
inertial range increases, the signed measure is called sign singular. In this case, changes 
in sign occur everywhere and in any scale considered in the limit of infinite Reynolds 
number. On the other hand, for a smooth field k = 0. 



3. Code and numerical simulations 

The data we use for the analysis stems from DNS of the incompressible Navier-Stokes 
equation with constant mass density, 

— +v-vv=-VP + ^V 2 v + f, (3.1) 

V -v = 0. (3.2) 

where v is the velocity field, p is the pressure (divided by the mass density), v is the 
kinematic viscosity, and f represents an external force that drives the turbulence. 

For the analysis, we consider six numerical simulations at different Reynolds numbers 
and spatial resolutions, with different forcing functions that inject either energy or both 
energy and helicity. The simulations are listed in Table [l] and described in more detail 

). Spatial resolutions range from 256 3 to 1024 3 
grid points, with Reynolds numbers R e — UL/v (with U the rms velocity and L the 
integral scale) ranging form ss 675 to « 6200, and Taylor Reynolds numbers R e = UX/u 
(where A is the Taylor scale) ranging from sa 300 to « 1100. Two forcing functions were 
considered: Taylor-Green (TG) forcing, which injects no helicity in the flow (although 
helicity fluctuations around zero develop as a result of nonlinearitics in the Navier-Stokes 
equation), and Arn'old-Beltrami-Childress (ABC) forcing, which injects maximum (pos- 
itive) helicity in the flow. The forcings were applied respectively in the shells in Fourier 
space with wavenumber kp = 2 and 3. All runs were continued for over ten turnover 
times after reaching the turbulent steady state. As a measure of the helicity content of 
each flow, we computed a relative helicity pn = H /{UfE) (where E is the total energy) 
averaged in time over the turbulent steady state of each run (see Table IT). 



in a previous paper (Mininni et al. |[2006 
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4. Results and analysis 

To study scaling laws of helicity fluctuations in the inertial range of helical and non- 
helical flows, we performed a signed measure analysis and computed the cancellation 
exponent n for all runs. Following Eq. (2.3 1, this can be done by computing %(7) and 
fitting x(l) ~ I ~ K in the inertial range. 

Although a detailed analysis of the simulations can be found in Mininni et al. ( 2006 1 , 
in Fig. [4] we show as a reference the energy spectrum of the two simulations at the 
larger Reynolds numbers attained (runs A3 and T3) . The spectra present a short inertial 
range followed by a bottleneck, as is usual in simulations of hydrodynamic turbulence. 
The runs also have a range with approximately constant energy flux, as well as a range 
of scales consistent with Kolmogorov 4/5 law when third-order structure functions are 
computed. As usual, these ranges have slightly different extensions depending on the 
quantity studied, and it will be shown that the cancellation exponent shows a scaling 
range comprised betwen I w 0.08 and w 0.7, which correspond roughly to wavenumbers 
k w 8 and 80, which is wider than but includes the inertial range of both runs. At the 
lower resolution, the inertial range in the energy spectrum and flux, or in the third-order 
structure function, is shorter and harder to identify. 

Fig. [2] (left) shows the signed measure of helicity as a function of the scale for runs Tl, 
T2, and T3 from top to bottom (non-helical, with increasing Reynolds numbers). As the 
Reynolds number increases, a wider scaling range can be identified. The values obtained 
for the cancellation exponent are k = 0.92 ±0.09, 0.75 ±0.03, and 0.73 ±0.01 respectively 
for runs Tl, T2, and T3. The values obtained and the dependence with Reynolds number 
support the idea that in the limit of vanishing viscosity, helicity fluctuations are sign 
singular, i.e., that fast changes in sign of helicity occur everywhere in arbitrarily small 
scales in the flow. This result is consistent with helicity cascading towards small scales, 



and with observations of the helicity distribution being highly intermittent (Chen et al. 
2003 b\ (albeit the previous studies are only for flows with net helicity) 



As mentioned earlier, the flows with ABC forcing are helical. For such flows the direct 
cascade of helicity and its intermittency has been studied before using spectral quantities. 
Although helicity fluctuations develop in these flows, given the predominant sign of helic- 
ity, fluctuations occur around the mean value and less changes in sign take place. Indeed, 
if the signed measure of helicity and the cancellation exponent is computed for run A3, 
we obtain k = 0.26 ± 0.04. This smaller value of the exponent is consistent with the less 
cancellations, and the larger error is associated to the fact that larger fluctuations are 
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Figure 2. Signed measure of helicity as a function of the scale for runs with increasing resolution, 
from top to bottom: TG forcing (left), and ABC forcing (right) after subtracting the mean value 
of helicity. Slopes with the cancellation exponent are given as a reference. 



observed in the signed measure as less events whith changes of sign are available. How- 
ever, we are interested in helicity fluctuations around the mean value, and to correctly 
consider these fluctuations the mean value of helicity H was subtracted from the helicity 
density in runs Al, A2, and A3 prior to the analysis. The resulting signed measures are 
also shown in Fig. [2] (right). In this case we also observe scaling of x(l) expanding over 
a wider range of scales as the Reynolds number increases, indicating that also for helical 
flows helicity fluctuations are sign singular. 

The results obtained for the helicity cancellations in helical and non-helical runs are 
similar. The cancellation exponents are K = 0.89 ± 0.07, 0.80 ± 0.03, and 0.72 ± 0.01 
respectively for runs Al, A2, and A3. At the highest resolution three ranges can be 
identified in both the A3 and T3 runs. At the smaller scales, viscous effects dominate 
and the flow is smooth. This results in less changes of sign in helicity and a shallower 
distribution of the signed measure. The slope tends to zero, as expected for a smooth 
flow, and the signed measure tends to one for the smallest resolved scale. At intermediate 
scales an inertial range is observed, and at larger scales the results are affected by the 
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1000 -200 -100 100 200 

R» H 

Figure 3. (a) Cancellation exponent k as a function of the Taylor based Reynolds numbers for 
the helical (triangles) and the non-helical (diamonds), (b) Histogram of helicity fluctuations for 
the T3 (solid) and A3 (dashed) runs. 



external forcing. In this scales, large fluctuations in the value of x(0 are observed as a 
result. As mentioned before, the inertial range of x(Z) is wider than the inertial range 
observed in the energy spectrum. This is consistent with observations of the helicity 
spectrum decaying slower than the energy spectrum near the dissipative range in helical 
turbulence ( Mininni et al. |[2006 l. 

The cancellation exponent obtained in all runs as a function of the Taylor based 
Reynolds number is shown in Fig. [3] (a). For both the helical and non-helical runs n 
decreases with the Reynolds number, and for the largest Reynolds numbers studied the 
value of k seems to saturate and is the same within error bars for both helical and 
non-helical flows. Note that because of the different forcing functions, the helical and 
non-helical runs have slightly different Reynolds numbers even when the same grid res- 
olution is used. The similarities between the scaling laws of the helical and non-helical 
flows suggests that statistics of the helicity fluctuations are the same in both cases. This 
is further confirmed by a histogram of the helicity fluctuations around its mean value for 
the T3 and A3 runs. As Fig. [3] (b) shows, turbulent fluctuations of helicity are similar, 
with small differences on the tails. 

Helicity fluctuations are sign singular. The fast oscillations in sign point to a highly 
intermittent quantity that fluctuates rapidly in localized structures at arbitrarily small 
scales for infinite Reynolds numbers. As mentioned in the introduction, the cancellation 
exponent can be related to the geometry of such structures, and under some conditions 



these relations can be rigorously derived (Vainshtein et al. 1994 1. The relations can be 



found also using simple geometrical arguments. In the following, we extend the analysis of 



Sorriso-Valvo et al. (2002) for the current density in magnetohydrodynamics to consider 



the case of kinetic helicity in hydrodynamic turbulence. We assume the helicity is spatially 
correlated (although not necessarily locally smoth) in D dimensions and uncorrelated in 
3 — D dimensions, where 3 follows from the dimension of the space. With this choice, 
a correlated helicity distribution has D = 3, and a completely uncorrelated helicity 
distribution has D = 0. For intermediate cases, we can estimate the signed measure of 
helicity as follows 



x(0 = E 

Qi(0 



/ d 3 xH(x) / I d 3 x\H(x) 
JQi(l) I JQ(L) 

[ d 3 xH(x) 

Jam 



(4.1) 
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Figure 4. Contour levels of helicity density in a slice of the T3 run. Only 1/4 of the slice is 
shown. Left: contours with i/(x) = 0, i.e., places where cancellations take place. Right: contours 
with H(x) = 12. Note how contour levels of zero helicity fill the space, while white patches are 
observed in the contour levels of non-zero helicity, as well as filamentary structures. 



where we used homogeneity to replace the sum over all subsets Qi(l) by (L/l) 3 (the total 
number of terms in the sum) times the integral over a generic box Q(l) of size I. Since 
the absolute values in the denominator prevent any cancellation, the integral below was 
approximated by the characteristic value L 3 H rms where H rms — (H 2 ) 1 / 2 is the rms value 
(jf the helicity density. 

In the inertial range, and still under homogeneity asumptions, the helicity follows some 
scaling law which can be represented by helicity structure functions (6H(s)) ~ s , where 
s = |s| represents a lineal increment. We can then split the integral over the domain Q(l) 
in integrals over subdomains with volume A 3 such that they separate the contributions 
from domains where the helicity is correlated of domains where the helicity is uncorrelated 



in space. Then, the remaining integral in Eq. (4.1 1 can be estimated as ( Sorriso-Valvo et 
aT|[2002 l 



I d 3 xH(x) -(H 2 ) 1 ' 2 / d D s r~) / d 3 - D s. (4.2) 

Q(0 JQ(l) VA/ JQ(l) 

The uncorrelated dimensions give a value proportional to the square root of their volume 
(1/ A)( 3-D '/ 2 . The correlated dimensions give a contribution proportional to (l/X) h+D , 
wich for h — (locally smoth helicity) is proporcional to their volume. Replacing this 



results in Eq. (4.1 1, we get that the signed measure of helicity x(J) scales as 



-h) 



(4.3) 



Thus the cancellation exponent n, the fractal dimension D, and the scaling exponent 
h are related through 

3-D 

K=—^--h. (4.4) 
A completely smooth field (D = 3) has k = h = 1, in agreement with the definition of 



k and of the first order (Holder) scaling exponent h for the field (see Bruno & Carbone 



(1997)). In a turbulent flow, the global helicity is observed to follow Kolmogorov scaling, 

1/3 ( |Kurien ||2003| |Mininni et aT]|2006[ ) . Then from 



and therefore we can assume h 

the value of k — 0.73 ± 0.01 it follows that helical structures are filamentary with fractal 
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dimension D = 0.873 ± 0.005, in good agreement with previous observations of helical 
vortex filam ents ( |Moffatt fc Tsinober |1992{ |Farge et a?.|2001| |Chen et aZ.|20036| |Mininni 
et al. [2006 ). As an illustration, Fig|4| shows a slice of the helicity density in the T3 run 



Filamentary structures can be identified, specially for the contour levels corresponding 
to large concentration of helicity. Structures less elongated were confirmed to correspond 
to filaments that cross the slice perpendicularly. 



5. Conclusions 

Signed measures of helicity were calculated for turbulent flows with different Reynolds 
numbers and helicity contents. For both helical and non-helical flows cancellation expo- 
nents were calculated, obtaining a value compatible with k = 0.73 ± 0.01 for helicity 
fluctuations for both cases at the largest Reynolds number considered. The scaling range 
increases with the Reynolds number, and the value of k obtained indicates helicity is sign 
singular in the limit of infinite Reynolds number. Simple geometrical arguments were used 
to link the cancellation exponent k to the fractal dimension D of the structures. To do 
this we assumed a single scaling exponent h for the lineal helical increments. Assuming 
h = 1/3, a value consistent with numerical results, we obtain that helicity fluctuations 
are filamentary with fractal dimension of D = 0.873 ± 0.005, also in agreement with 
observations. 

The authors acknowledge support from Grants No. UBACYT X468/08 and PICT- 
2007-02211. PDM acknowledges support from the Carrera del Investigador Cientffico of 
CONICET. 
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